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In this work, a new functional is introduced to treat pairing correlations in finite 
many-body systems. Guided by the projected BCS framework, the energy is written as 
a functional of occupation numbers. It is shown to generalize the BCS approach and 



r^ i to provide an alternative to Variation After Projection framework. Illustrations of the 

^ ' new approach are given for the pairing Hamiltonian for various particle numbers and 

coupling strengths. In all case, a very good agreement with the exact solution is found. 
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f~^ ■ 1. Introduction 

lO ' The Nuclear Energy Density Functional (EDF) theory is anticipated to provide a 



unified framework to describe nuclear structure and reactions. Important efforts are 
now being made to improve the predicting power of EDF. An important aspect of 
current EDF is the possibility to break some of the symmetries of the nuclear Many- 
Body problem in order to grasp specific correlations with rather simple functionals. 
p% ' A typical example is provided by pairing correlations that are approximately treated 

jrt ' by relaxing the particle number conservation, like in the BCS or HFB theory. Such 

theories are expected to be rather effective in the large particle number limit but 
miss important effects as the particle number and/or coupling decreases. Projection 
techniques that restore the particle number are expected to significantly improve 
the description of pairing in that casea^. Figure [T] illustrates the energy for the 
"picket fence" pairing Hamiltonian 1^, 

obtained in the BCS, and projected BCS [PBCS] frameworks when the variation is 
made before (Variation After Projection [VAP]) or after the variation (Projection 
After Variation [PAV] ) . The projection clearly improves the description of correla- 
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tion and is almost indistinguishable from the exact result when VAP is performed. 
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Fig. 1. Difference between the exact total energy and the Hartree-Fock energy (solid line) ob- 
tained in the picket fence Hamiltonian with constant level spacing Ae for A = 8 and yl = 16 
particles. In both case, the BCS (dashed), PAV (open triangle) and VAP (open circles) results are 
displayed. 



Projection on particle number and/or angular momentum is becoming now a 
standard tool of nuclear EDF and extensive calculations along the whole nuclear 
chart is within reach. Recent discussions have however pointed out tha t projection 
should be handled with care when combined with functional theory^l^. In particu- 
lar, specific correction to th e so-called self-interaction problem in "GCM like" cal- 
culations are necessar y I I I In addition, even with proper corrections, projections 
especially made before the variatiorM^ remain very heavy numerically and their use 
for other problems than nuclear structure, like for instance the dynamics or ther- 
modynamics of nuclear systems, seems difficult. The aim of the present work is to 
propose a new functional for pairing able to directly grasp correlations generally 
incorporated through projection by specific dependence on the the natural occupa- 
tion numbers Ui and orbitals {vi}, i-e. p = X]i |'/'i)"-i(Vi| where p denotes the exact 
one body density matrix. 



2. Density Matrix Functional for pairing 

Our starting point is the projected BCS state given by: 
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(2) 



where {a], at} denotes pairs of time-reversed states. It could easily be shown 
that this state can be obtained by projection of the BCS state, |BCS) — J|(l + 
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XiaJaJ) |0)|—), onto good particle number A — 2N. Note that here N denotes the 
number of pairs. The PBCS energy for the Hamiltonian ^ reads 

{N\H\N) _ 



{N\N) 



Sijii, dj) = 2 ^ EiUi - .9 ^ Cij, (3) 



where n,; and Cij are the occupation probabihties and correlation matrix elements 
defined through 

_ {N\a\a.\N) {N\h\b,\N) 

'~ {N\N) ' '■'■" {N\N) ■ ^> 

According to the definition ^, rii and Cij and therefore £{ni, Cij) can eventually be 
written as an explicit functional of the {xi} parameters, which turns out to be too 
complicated for a direct use as variational parameters. After tedious manipulation, 
it has been shown in ref. '^, that these parameters can approximately be written as 
a functional of the {rii} through; 



X, 



|2 



(ao ~ ain.i) (5) 



with 



and 



ar = l(l + .. + .^ + ... + 4-)4l3| (6) 



ao = l + ^^^^(l + 2s2 + --- + (iV-l).r2) = l + (52-.3)^, (7) 

and Sp — l/NY^ nf . With these expressions, we can now write the correlation as a 
functional of single-particle occupancies: 

^/(ao - aini){ao - aiUj) 
{ao - ai{ni + Uj - n^rij)} 



C,j = \/n^{l - n,)nj(l - rij) x -f^ ^^^ ' TTT^T- i^) 



The present functional, based on the PBCS trial state, corresponds to a general- 
ization of the BCS ansatz that approximately account for particle number conser- 
vation. The BCS limit is obtained for qq = 1 and oi = leading to the well know 
expression Cij — y^ni(l — ni)nj{l — rij). It is worth to mention that a completely 
different expression is obtained in the weak coupling limit (Hartree-Fock limit) for 
which S2 = S3 = 1 and Cij — )■ ^Juijij. 

Illustrations of the new functional accuracy are given in figure [5] where the 
energy is obtained by a direct minimization of ([3|) using expression ([8]) . A very good 
agreement is obtained for any particle number and coupling strength. 
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Fig. 2. Same quantity as in Figure [T] displayed as a function of the coupling strength for the 
exact (solid line), BCS (dashed line) and the new functional (filled circles) case. 

3. Summary 

Following the Density Matrix Functional Theory spirit ^^^ the possibility to use 
functionals of natura l orbitals and occupancies has been recently introduced in 
nuclear physics^ I i In the present work, a new functional for pairing is introduced 
that correct the BCS theory for particle number conservation effects. The present 
framework provide an alternative way to perform VAP calculations and is expected 
(i) to greatly simplify such calculation compared to direct use of projection operator, 
(ii) to directly give access to physical quantities like natural occupancies, (iii) to be 
particularly suitable for EDF based theories. 
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